Lagrangian tracer dynamics in a closed cylindrical turbulent convection cell 



O 

(N 

C 
3 



Mohammad S. Emran 1 and Jorg Schumacher* 1 

1 Institut fur Thermo- und Flutddynamik, Technische Universitat Ilmenau, Postfach 100565, D-98684 Ilmenau, Germany 

(Dated: June 7, 2010) 

Turbulent Rayleigh-Benard convection in a closed cylindrical cell is studied in the Lagrangian 
frame of reference with the help of three-dimensional direct numerical simulations. The aspect 
ratio of the cell T is varied between 1 and 12, and the Rayleigh number Ra between 10 7 and 
10 9 . The Prandtl number Pr is fixed at 0.7. It is found that both the pair dispersion of the 
Lagrangian tracer particles and the statistics of the acceleration components measured along the 
particle trajectories depend on the aspect ratio for a fixed Rayleigh number for the parameter range 
covered in our studies. This suggests that large-scale circulations present in the convection cell 
affect the Lagrangian dynamics. Our findings are in qualitative agreement with existing Lagrangian 
laboratory experiments on turbulent convection. 
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I. INTRODUCTION 

The motion of a fluid can be described either in the Eu- 
lerian or Lagrangian frame of reference. In case of fluid 
turbulence, the Lagrangian point of view brought inter- 
esting and new insights on the small-scale structure and 
statistics (see e.g. Yeung [l|, Toschi and Bodenschatz 
Q for comprehensive reviews). These studies included 
intermittency of temporal velocity increments and accel- 
erations [l-Q or the geometry of particle tracks @ . Fur- 
thermore, two-particle and multi-particle dispersion have 
been analysed and compared with classical predictions by 
Batchelor Q and Richardson 

Almost all of the experimental and numerical investi- 
gations in turbulent thermal convection have been con- 
ducted in the Eulerian frame. This includes studies on 
the turbulent heat transfer and large-scale circulations 
d, [l(| as well as on small-scale structures and dynam- 
ics of thermal plumes [TlT - fll| . Recently, first Lagrangian 
laboratory experiments on convection were conducted in 
a closed cylindrical vessel [l4| and numerical simulations 
were carried out in Cartesian slabs with periodic side 
walls and free-slip boundary conditions at the top and 
bottom [l5|, [lj| . The focus of these studies was on the 
local variations of the heat transfer which can be mea- 
sured along the trajectories of the Lagrangian tracers. In 
the numerical studies, these results could be connected to 
the local tracer accelerations. In the laterally infinitely 
extended layer (i.e. the configuration with periodic side 
walls), slightly less intermittent vertical accelerations a z 
were detected in comparison to the lateral ones, a x and 
a y . This was manifested in the sparser tails of the prob- 
ability density function (PDF) of a z . 

In the present work, we want to take these numerical 
Lagrangian studies to the next level. We report investi- 
gations in a closed cylindrical cell with no-slip boundary 
conditions at all walls. This is the characteristic setup 



for almost all laboratory experiments. It will allow for a 
closer comparison with the work by Gasteuil et al. 14 1. 
In addition, we want to explore systematically how the 
finite size of the cell affects the Lagrangian dynamics, 
such as the long-time behaviour of the two-particle dis- 
persion. The latter point is also motivated by recent Eu- 
lerian studies, in which a systematic dependence of the 
turbulent heat transfer on the aspect ratio T is found for 
small and moderate values of V [IT} • The dependence is 
in line with morphological changes in the large-scale cir- 
culation (LSC), e.g. from one-roll to multi-roll patterns, 
in the convection cell. Our boundary conditions will dif- 
fer from those in Refs. jl{iLll6| . This difference might also 
affect the statistics of the acceleration components. It is 
thus another open question if a different velocity bound- 
ary layer structure has an impact on the intermittency 
of lateral and vertical acceleration components. 

The outline of the paper is as follows. In section II we 
present the equations of motion, the numerical model and 
the tracer advection scheme for the non- uniform cylindri- 
cal mesh. Section III is focussed on the particle dynamics 
and local heat transfer followed by studies of the pair dis- 
persion in Section IV and the aspect ratio dependence of 
the acceleration statistics in Section V. Finally, we sum- 
marize our findings and give a brief outlook. 



II. EQUATIONS OF MOTION AND 
NUMERICAL MODEL 



A. Boussinesq equations 

We solve the Boussinesq equations for Rayleigh- 
Benard convection numerically with a second-order fi- 
nite difference scheme in cylindrical coordinates {4>, r, z) 
[l9l [2fjj . The equations are given by 



du 

~dt 
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(u-V)u = -Vp + vV 2 u + agTe z , (1) 
V-u = 0, (2) 
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193 x 97 x 128 


10 7 
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2.26 
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257 x 165 x 128 


10 7 
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1.43 


14 
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361 x 257 x 128 


10 7 
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1.32 


14 
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401 x 311 x 128 


10' 
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1.06 


14 
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601 x 401 x 128 


10 7 


12 


1.01 


14 


6 


271 x 151 x 256 


10 s 


1 


1.73 


20 


7 


361 x 181 x 310 


10 9 


1 


1.03 


17 



TABLE I: Simulation parameters for various Ra and Y with 
a fixed Pr = 0.7. The grid resolution, the Rayleigh number 
Ra, the aspect ratio Y, the ratio nrj/A, and Nbl, the number 
of horizontal grid planes inside the thermal boundary layer 
thickness, are given. Nbl satisfies the criteria discussed in 
[l8| very well. 




FIG. 1: Sketch of a 3D grid cell in the cylindrical coordinates, 
in which the positions of the velocity components, the tem- 
perature and the pressure are shown on the staggered grid. 



dT 
~dt 



+ (u-V)T = kV 2 T, 



(3) 



where p(x,t) is the (kinematic) pressure, u(x,t) the ve- 
locity field, a the thermal expansion coefficient and g 
the gravity acceleration. The velocity field has a no-slip 
boundary condition at all walls. The temperature field 
is isothermal at the top and bottom plates and adiabatic 
at the side wall. Following Refs. [lj, [l7|, we have cho- 
sen the DNS grid such that the global maximum of the 
geometric mean of the grid spacings in ip, r and z, i.e. 



A = max[^/rA A r A 2 



(4) 



cylindrical coordinate frame. The spatial discretization 
is performed on a staggered grid and solved by a second- 
order finite difference scheme [III [2(| . The pressure field 
p is determined by a two-dimensional Poisson solver [22[ 
after applying one-dimensional fast Fourier transforma- 
tions in the azimuthal direction. The time advancement 
is done by a third-order Runge-Kutta scheme. The grid 
spacings are non-equidistant in the radial and axial di- 
rections. In the vertical direction, they correspond to 
Tschebycheff collocation points. Some parameters and 
the corresponding grid resolutions are given in Table |TJ 



satisfies the resolution threshold by Grdtzbach [2l| . 
namely nr/fA > 1 where rj is the Kolmogorov dissipa- 
tion length which is calculated via the mean energy dis- 
sipation rate (e) = jpc(Nu — l)RaPr~ 2 . The details are 
provided in Tab.l. 

The resolutions of the runs are exactly the same as in 
[l7| . Note that the turbulence is strongly inhomogeneous 
with respect to the vertical direction z. This results in 
amplitudes of the energy dissipation rate which are by 
orders of magnitude larger close to the plates than in the 
bulk. A local Kolmogorov scale is therefore much larger 
in the bulk as in the boundary layers. It motivated us 
in [TtJ to relate the grid spacing to i](z) = v z / A / (e(z))A,t 
rather than r\. 

The dimensionless parameters and their range of val- 
ues in our numerical study are: the Rayleigh number 
Ra = ag(Tb ottom — T top )H 3 /(vk) between 10 7 and 10 9 , 
the Prandtl number Pr = v/n = 0.7 and the aspect ra- 
tios r = D/H between 1 and 12. D is the cell diameter 
and H is the height of the cell. 



B. Numerical integration scheme 

We conduct direct numerical simulations (DNS) of the 
Boussinesq equations. The equations are solved in a 



by 



Tracer advection 



The equations of Lagrangian tracer motion are given 



dx p (t) 
dt 



u(x p ,t) . 



(5) 



Here u = {u c f,{4>p),u r (r p ),u z (z p )) T is the velocity of the 
particle at position x p — {<j} v {t),r p {t),z p {t)) T with the 
symbols <f>, r and z corresponding to the azimuthal, ra- 
dial and axial directions, respectively. As discussed 
before, the Eulerian fields are well resolved on a stag- 
gered and non-uniform grid in the simulations. The ve- 
locity components, as shown in Fig. [TJ are then given on 
particular faces of the grid cells only, while the Eulerian 
temperature field is always centered in the grid cell. All 
these fields are required at the Lagrangian particle po- 
sition to advance (J5J) and study local heat transfer. We 
therefore applied an interpolation scheme for the veloc- 
ity components which makes direct use of the staggered 
setup. The interpolation is then of second order and the 
maximal error can be estimated to be of the order of 
C(A 2 ) (see (@|). Given the Eulerian velocities at the 
mid points of the edges of the mesh cell (see Fig. |5J , we 
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u r (i,j-i/2,k) >e 



u z (i,j,k+1/2) 

X 



gives 



W r j.1/2 u r( r p) w r j+1/2 



Ar, 



-X- 



U r ,L - r L - 



Au r 
Ar 



(9) 



Now integrating both sides of ^ with respect to t gives 
the r-position of the particle as 



X u r (i,j+1/2,k) 



r p (t) = rx, - w r ,_ 



> r u z (i,j,k-1/2) 

FIG. 2: Position of the velocity components on a 2D staggered 
grid var — z plane, which is the mid section of a typical 3D cell 
as shown in the top panel of the figure. The indices i, j and 
k correspond to the azimuthal, radial and axial directions, 
respectively. The particle velocity u r (r p ) at position r p is 
obtained by the linear interpolation of the Eulerian velocities 
u r (i,j — \,k) and u r (i,j + \,k) with the widths w r ._i and 



w r , i as given in (|6} . 

3 + 2 



J-2 



calculate the corresponding Lagrangian components at 
an arbitrary position inside a cell as 



u r (r p ) 

u<f>(<P P ) 
u z (z p ) 



w r _±u r (i,j + |,fc) + w r j+ iu r (i,j - |, k) 
A r (i, j, fc) 

w^iU^i + \,j,k) +wf + iu 4> (i - fc) 

w z k _iu z {i,j,k+ i) + w z k+ iu z (i,j,k- §) 
A 2 (i,j,fc) 



(6) 



Indices fc denote the grid positions with respect to 
cj>,r,z, which run from 1 to A^, AT r and A^, respectively. 

In order to progress further, we rewrite the three or- 
dinary differential equations (we provide here details of 
the equation in r-direction only) [23| as 



dr p 
~dt 



u r (r p ) 



Tr - r L 



tl 



tr - r L 



Ur,L ■ 



where M r ,i? and it r ,L are the Eulerian velocities at face 
positions tr and rj,, respectively and r p is the particle's 
position in between. The subscripts L and R refer to the 
left and right faces of a 3D-cell (see Fig. [IJ , respectively. 
Rearranging ([7]) yields [23j 



dr p 
dt 



Ai 



p Ar 



= U r ,L - T L ■ 



Au r 
Ar 



(8) 



with Au r = u Tt R — it r ,z, and Ar = rR — rL- Multiplying 



both sides of © with the integrating factor exp(- 



Am, 



t) 



At t = t 



r P (t ) = r L - u r j 



Ar 

Au r 



Ar 
Au r 



C> 



C> 



and thus 

Sr = r p (t) - r p (t ) = C \e^ t - 

From Eq. (|lip . with r p (to) = r , and Eq. 
the constant C follows to 

Ar A« r . 
C = u r . — e -^*° . 
Au r 

and thus the radial increment for 8t = t — to is 
A 



(10) 

(11) 

• (12) 
at t = to 



St = u r ,o-^- (e^ st - l) 

AUr V / 



(13) 



(14) 



Expanding the exponential term, this scheme is first or- 
der in time for the leading contribution. Similarly, for the 
azimuthal and axial directions, the distances are given by 



ro Au^ 



oz = u z .o- I e 



'At 



(e^fr* - l) 



(15) 
(16) 



The equations of the particle trajectories in the La- 
grangian frame are 



r p (t) = r p (t Q ) + Sr, 
4> P (t) = 4> P {to) + Sep, 
z p (t) = z p (t ) + 5z . 



(17) 
(18) 
(19) 



The temperature field T is stored at the cell center 
of the staggered grid as shown in Fig. [TJ At first, we 
interpolate T at the vertices of the Eulerian grid and then 
apply trilinear interpolation to calculate the Lagrangian 
temperature at a known particle position. 



III. TRACER DYNAMICS IN THE CELL 

A. Tracer trajectories and large-scale circulation 

For the Lagrangian simulations, we have selected cylin- 
drical cells which are given in Tab. U keeping the Prandtl 
number constant (Pr = 0.7). We advect N p ~ 10 5 tracer 
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Run t/tf 


Nu E 4= a 


Nu L ± a 


N p 


1 


150 


16.734=0.08 


16.774=2.95 


104160 


2 


150 


16.06±0.05 


15.864=1.37 


81024 


3 


150 


16.364=0.03 


15.70±0.57 


202860 


4 


150 


17.44±0.02 


16.00±0.54 


208425 


5 


150 


17.49±0.03 


16.064=0.54 


272811 


6 


150 


32.21±0.32 


31.754=4.29 


112320 


7 


139 


64.31±0.64 


65.124=19.53 


179850 



TABLE II: The total integration time t/tf, the Eulerian and 
Lagrangian values of the Nusselt number plus errorbars Nue 
and Nul, and the total number of tracers N p are listed. 



particles in the cell (see Tab. |II]for more details). They 
are initially seeded in the whole cell and integrated in 
time for more than 150 time units t/tf where tf = H/U f 
is the free fall time. Here, Uf = yj ga(Tbottom — T top )H. 
The statistical analysis is thus conducted over 3000 in- 
dependent Lagrangian samples, which are separated by 
At = 0.05£/ which corresponds to At 

10 7 , 10 8 and 10 9 , respectively. 



0.25r n , 0.33t,, 



and 0.43r^ for Ra 




FIG. 3: (Color online) Three-dimensional trajectory of a La- 
grangian tracer particle for the simulation at Ra = 10 7 in the 
cell with aspect ratio T = 1 . The tracer track is displayed 
for 100 free-fall time units tf. 



We visualize the time-evolution of the trajectory of a 
Lagrangian particle in Fig.[3]for a simulation at Ra = 10 7 
in the Rayleigh-Benard cell with aspect ratio T = 1. It 
seems that this particle travels preferentially along the 
large-scale circulation (LSC). Fig. |4] plots the time se- 
ries of the vertical position and various quantities along 
the same particle path as shown in Fig. [3J The results 
are qualitatively in agreement with the Eulerian mea- 
surements [24[ and the Lagrangian measurements [lij ]. 
Sharp spikes of the temperature signal appear due to ris- 
ing and falling plumes close to the hot and cold plates of 
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FIG. 4: (Color online) Time series of various quantities along 
the trajectory of a single Lagrangian tracer particle: (a) the 
vertical position z(t) normalized by the cell height H , (b) the 
vertical velocity u z (t), (c) the temperature T(t), products of 
velocity components and temperature (d) u z T, (e) u y T and 
(f) u x T. The dotted lines in some of the panels show the 
corresponding value of the time average which enters the La- 
grangian definition of the Nusselt number as given in Eq. (I21|l . 
Data are taken along the tracer track which is shown in Fig. 

El 



the cell. They are correlated with the maxima and min- 
ima in the vertical position. The vertical velocity u z (t) 
has zero mean. The convective heat fluxes are shown in 
Fig. H)Jd)-(f). The signals are highly irregular with sharp 
intermittent peaks. The factor %/ RaPr is multiplied to 
all flux terms u x T, u y T and u z T in order to get the values 
of the dimcnsionless local convective heat flux. Similar to 
Gasteuil et al. [l4| we observe strong fluctuations about 



5 



the mean value which are 4.8, — 1.7 and 23.4 for u x T, 
u y T and u z T, respectively. It demonstrates that the ver- 
tical flux u z T is dominantly responsible for the net heat 
transfer in the cell. 

The time series of the vertical tracer position z(t) 
shows a rather regular oscillatory form which is caused 
by the mean wind in the cell. One can therefore extract 
a typical loop time from all the time series. Motivated 
by the analysis of Sreenivasan et al. 25] for mean wind 
reversals, we have therefore performed a statistics of the 
zero-crossings for £,•(*) with i = l..N p . The mean time 
interval between two subsequent zero crossings for each 
tracer track has been determined, and a subsequent La- 
grangian ensemble mean, (•}/, has been taken. The lat- 
ter, which we denote as (T zc )l, gives a characteristic loop 
time of the tracers, X/ = 2(T zc )i J . For a fixed Rayleigh 
number Ra = 10 7 , Tj = 19.4*/, 17.8*/ and 19.0*/ for 
r = 1, 3 and 8, respectively. 

Translating into a convective time unit t c = H/(u 2 )y^ 
(seeRef. Q3 for a discussion) gives TJ ~ 3.7* c , 3.4* c and 
3.6 * c for r = 1,3 and 8, respectively. For a fixed as- 
pect ratio r = 1, we find T; = 19.6*/ at Ra = 10 s and 
T; = 20.6*/ at Ra = 10 9 . We see that the characteristic 
loop time of the tracer is slightly varying but no clear 
trend when T is varied and Ra is fixed. It seems, there- 
fore, to be independent of the presence of a single-roll or 
multi-roll LSC. In case of the fixed aspect ratio, a slight 
increase of Tj with increasing Ra is detected. We recall 
that the velocity fluctuations amplitude decreases with 
increasing Ra as reported for example by Verzicco and 
Camussi [20]. Our observation of a slightly longer loop 
time seems to be consistent with a reduced amplitude of 
velocity fluctuations with increasing Ra. On the basis of 
the present data, it is however difficult to draw a firm 
conclusion about the robustness of this trend. 



B. Heat transfer in the Lagrangian frame 

Similar to the Eulerian case, where the Nusselt number 
is given by 

Nu E = l + ^-{u z T) v>t , (20) 

with (-)v,t being a volume and time average, the Nusselt 
number in the Lagrangian frame, Nul, is given by [l6[ 

Nu L = \ + -?—{u z T) L<t = {Nu L {t)) t , (21) 

where the symbol {-)L,t denotes an averaging over all 
trajectories and time. We compute the time series of the 
heat transfer in the Lagrangian frame and compare with 
its value in the Eulerian frame in Fig. [5] for a simulation 
at Ra = 10 7 in the cell with aspect ratio T = 1. The 
time averaging of Nu L {t) = 1 + H/(nAT)(u z (t)T(t)) L 
curve from * = 40 to 200 yields Nul = 16.77, which 
almost exactly matches the corresponding value of the 
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FIG. 5: (Color online) Graph of NuL(i) which is averaged 
over ~ 10 5 tracer paths. The temporal mean value of Nul = 
{Nu L (t))t (see Eq. IpT])) between time t = 40-200 is Nu L = 
16.77, which almost exactly matches with the corresponding 
Eulerian value Nue = 16.73 as listed in Table [TT] 



Eulerian simulation, which was Nue — 16.73 as pro- 
vided in Tab. [TTJ Similar result is achieved for the simu- 
lation conducted at Ra — 10 8 in the same cell, for which 
Nu = 32.21 and Nul = 31.75 with a deviation of 1.5% 
from Nue- The results also validate our interpolation 
schemes and the time step sizes for the Lagrangian sim- 
ulations, which are accurate enough to reproduce flow 
quantities of the corresponding Eulerian simulations in 
turbulent convection. In contrast, recent experiments by 
[HI obtained a value of Nul almost twice as large as the 
Eulerian case. This discrepancy is due to the fact that 
their [l4| mobile sensor traversed preferentially along the 
mean flow circulation path and hence the contribution 
from the rest of the volume was missing in their Nusselt 
number measurement. 

At the beginning of the NuL(t) curve in Fig. O there 
is a large overshoot from mean, which generally appears 
when the majority of the particles follow the LSC path 
or the particles are not seeded uniformly in the flow. We 
also observe strong oscillations in NuL{t) for all the sim- 
ulations listed in Tab. [TT] Strong fluctuations around the 
mean result in large standard deviation a. The value of 
a decreases with increasing T. When the number of par- 
ticles was increased (case: N p 2.5 x 10 5 , Ra = 10 7 and 
r = 1; the results are not shown here), the convergence 
of NuL{t) was achieved within a shorter period of time. 
It can be concluded that the number of tracer particles 
should be sufficiently large and they should be uniformly 
seeded in the domain for faster statistical convergence of 
Nu L (t). 
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FIG. 6: Pair dispersion as in Eqns. (|22l) (I24|l : (a) aspect ratio 
r = 1, (b) r = 3, (c) pair dispersion normalized by {e)t 3 for 
three initial separations, r/, 2r\ and 4rj, with r\ the Kolmogorov 
length and (d) the autocorrelation functions of the lateral and 
vertical velocities for F = 1 as given in Eqns. (|25p and (|26[1 . 
respectively. Fhe Kolmogorov time t v (see Eq. (|29[) ) and the 
lateral Lagrangian time T { L L) (see Eq. pgjll are marked on 
the time axis in (a) and (b). Fhe horizontal dashed lines in 
(a) and (b) denote D 2 = F 2 H 2 . Fhe vertical dashed line in 
(d) represents the characteristic time required by a particle 
to travel the cell of height H. Fhe analysis is conducted at 
Ra = 10 7 with 4.16 x 10 s and 3.24 x 10 8 statistical samples 
for T = 1 and T = 3, respectively. 
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FIG. 8: (Color online) Fime evolution of the lateral pair dis- 
persion as a function of the aspect ratio V. Data are for 
Ra = 10 7 . Fhe lowermost curve is for F = 1 and the upper- 
most for T = 12. 



IV. LAGRANGIAN PAIR DISPERSION 

The pair dispersion measures the relative separation of 
a pair of tracer particles traversing along their trajecto- 
ries and is defined as 



R 2 (t) 



[x p {t] - x p (t )Y 



(22) 



i 0.5 
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— P2 






P3 
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— P5 
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FIG. 7: (Color online) Relation between the vertical coordi- 
nate and the anti-correlation of vertical velocity. In the top 
panel, we show the 2(t)-component of 6 particle trajectories, 
which were seeded close to the bottom plate. Fhe correspond- 
ing product u z (Q)u z (t) along each trajectory which enters 
into the definition of the Lagrangian integral time is shown in 
the bottom panel. 



Since convective turbulence is inhomogeneous, we decom- 
pose the distance vector into two parts - the lateral part 
Rxyit) and the vertical part R z {t). Their contributions 
to the total dispersion R 2 (t) are given by 



Rl v (t) = ([x p (t)-x p (t )] 2 + [y p (t)-y p (t )] 2 ) L (23) 



R 2 z {t) = ({z p (t)~z p (t )} 2 }i 



(24) 



This decomposition is chosen in order to relate our re- 
sults to the findings in [ID, [l6j]. We also compute the 
autocorrelation functions of the velocity components as 



C xv (t) 



U X y(t)) L ,t 



(\u x y\ 2 ) L ,t 
(u z (t + T)u z {t)) Ltt 



(u 2 z ) L ,t 



(25) 
(26) 



Here C xy (r) is the autocorrelation coefficient for the lat- 



eral tracer velocities u xy = u x e x 



u y e v , with e x and e y 



are the unit vectors in x and y directions, respectively, 
C z {t) is for the vertical velocity u z , and t is the time 
lag. The vertical and lateral Lagrangian times, T[ and 
T^ y are obtained by integrating C z (t) and C xy (r) as 



(27) 
(28) 



Tl = \ C z (r)dT 
Jo 

poo 

T X L V = / C xy (r)dr, 
Jo 
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T 1 3 5 8 
r„ 0.211 0.215 0.213 0.206 
Tl y 0.68 0.98 1.21 1.16 



TABLE III: Kolmogorov time scale and lateral Lagrangian 
integral times. Data are for Ra = 10 7 and different aspect 
ratios. 



and the Kolmogorov time t v is given by 




Figures |U[a) and (b) plot the total dispersion, R 2 , as well 
as the lateral, R 2 y , and vertical, R 2 , contributions. Two 
aspect ratios are selected, namely T — 1 and 3 for simu- 
lations at Ra — 10 7 . Initially, the quantities grow with 
R 2 ~ t 2 , which corresponds to the ballistic dispersion. 
After the ballistic growth, a transition to a regime with 
a growth law of R 2 ~ t 3 occurs. Similar observations 
were made for fluid turbulence I26ll and for convective 
turbulence in an extended layer [l5l [H[ . The range is 
however too short to conclude if this a Richardson-like 
regime 8] or not. The horizontal dashed lines mark the 
square of the cell diameter, the limit which cannot be 
exceeded by pair separation. The finite size of the cell 
suppresses also an eventual Taylor diffusion regime with 
R 2 ~ t and as a result the dispersion levels off, which 
is different from the Cartesian cells with periodic side- 
walls [TH, [l6{ . The Kolmogorov time T n and the lateral 
Lagrangian integral time T L V as a function of the aspect 
ratio are listed in Tab. IIIII While the Kolmogorov time 
remains nearly unchanged with an increase of T, the lat- 
eral Lagrangian integral time increases by almost 100% 
for an increase of F from 1 to 5. This would be in line 
with a stronger correlation of the lateral velocity compo- 
nents which could arise in the multi-roll LSC case. 

Recent numerical Lagrangian studies in homogeneous 
isotropic turbulence by Sawford et al. [29j j yielded a sen- 
sitive dependence of the pair dispersion on the initial 
pair separation. This was confirmed in the case of con- 
vective turbulence in an extended layer |16| . In addi- 
tion it was found there that the intermediate evolution 
of the dispersion in thermal convection depends on the 
seeding height of the tracer pairs since convective turbu- 
lence is inhomogeneous in the vertical direction. Here, 
we repeat these studies. In Fig. E]Jc), the pair dispersion 
normalized by (e)t 3 is plotted for three different initial 
separations, namely r], 2r] and 4r], with r\ = ^ 3 / 4 /(e) 1 / 4 
the Kolmogorov length. The figure shows that a plateau, 
i.e. R 2 ~ t 3 , is reached for the smallest initial separation 
only. 

The autocorrelation functions of the lateral and ver- 
tical velocities are shown in Fig. [6jd) for one run. The 
results are in agreement with [15| . The vertical veloc- 
ity is strongly anticorrelated to the lateral one. This 
strong anticorrelation is in line with the rapid upward 
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FIG. 9: (Color online) The probability density functions of 
the acceleration components for different Rayleigh numbers 
at r = 1. Amplitudes are normalized by the corresponding 
root-mean-square value, (a) Component a x . (b) Component 
a z . 



and downward motions as indicated by the signals of ve- 
locity components along a tracer path in Fig. |4l in which 
fluctuations in u z (t) are stronger than those in u x (t) and 
u y {t) (time signals of u x and u y are not shown here). In 
order to demonstrate this more quantitatively, we plot in 
Fig. [7] the vertical displacement of the Lagrangian par- 

(z) 

tides and the corresponding integrand of definition T£ 
(see Eg. |2"8)) ) for 6 tracers that are initially seeded close 
to bottom boundary layer. The pronounced minimum of 
w z (0)m z (t) coincides with the reversal point of z(t) close 
to the top boundary layer. 

Fig. [8] shows the dependence of the long-time behavior 
of the lateral dispersion on the aspect ratio. Only for 
the largest aspect ratio (r — 12) a diffusion limit with 
R\ y (t) ~ t is obtained. It underlines an important differ- 
ence in contrast to configurations with periodic boundary 
conditions in the lateral directions. In our case, the dis- 
persion is additionally constrained by the sidewalls for 
small- aspect-ratio systems. 
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FIG. 10: (Color online) The convergence of the probability 
density functions of one lateral and the vertical acceleration 
components for different Ra = 10 7 and Y = 1. Data are for 
2 x 10 8 , 4 x 10 s , 10 9 , and 2 x 10 9 sample events. 
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FIG. 11: (Color online) The probability density functions of 
the acceleration components, (a) Component a x . (b) Compo- 
nent o z . The uppermost graph is for the largest aspect ratio 
in both panels. 



V. ACCELERATION STATISTICS 

Finally, we report the probability density func- 
tion (PDF) of acceleration components along three- 
dimensional trajectories. In Fig. |9] we compare data for 
three different Rayleigh numbers with fixed aspect ratio 
r = 1. Both lateral and vertical acceleration compo- 
nents show no systematic trend with Ra. The stretched 
exponential form is similar to those reported in labora- 
tory experimen ts p -a and DNS simulations of turbulent 
convection flEl fl6|7 

A critical issue of numerical and experimental La- 
grangian studies - not only in turbulent convection - is 
the statistical convergence, which requires to determine, 
for example, a fourth-order moment, such as the flatness 
Fiflk)- In Lagrangian turbulence, the situation becomes 
particularly problematic since the tails of the PDFs are 
found to be strongly scattered [3(| [M[ . The statistical 
uncertainty is also manifested in the scattered tails of the 
PDFs in Fig. [9] although more than 2 x 10 8 data points 
have been included. In Fig. I1Q[ we ran the simulation 
for significantly longer period such that the number of 
samples became one order of magnitude larger. The tails 



become less scattered and more extended with increasing 
number of samples. The data indicate that the conver- 
gence in the vertical direction is slower than in the lateral 
direction. Nevertheless, the scatter for the fourth order 
moments remains so strong that the flatness is obtained 
with a big error bar only. The flatness is defined as 

F{a k ) = (30) 

\ a k)L.t 

with k — x,y and z. We found that F(a z ) = 89 ± 18 for 
4 x 10 9 samples. 

The dependence on the aspect ratio T is demonstrated 
in Fig. [TT] It can be seen that the tails of the PDF of 
a x grow with increasing aspect ratio, which implies that 
less constrained lateral motion causes larger acceleration 
magnitudes. Meandering of the tracers among different 
LSC rolls pjj will probably contribute to larger ampli- 
tudes in the acceleration as well. Interestingly, a similar 
but weaker trend holds for the vertical acceleration com- 
ponent as well. 
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VI. SUMMARY AND DISCUSSION 

We have presented DNS studies on the Lagrangian 
tracer particle dynamics in turbulent Rayleigh-Benard 
convection. The main objective of this work was to in- 
vestigate how the finiteness of the cell size in convection 
impacts the motion, dispersion and acceleration of the 
tracer particles. We developed an interpolation scheme 
for the staggered computational mesh, which was used 
in the Eulerian simulation. Our interpolation scheme re- 
produces accurately the global heat transfer in the La- 
grangian frame Nul, which matches the corresponding 
value in the Eulerian frame Nue- We detected however 
large fluctuations of Nul($) about the mean value since 
the majority of the tracers is moving along the mean cir- 
culation path in the cell. Our studies suggest that either 
a large number of particles or a very long time integration 
period are necessary for faster convergence of NuL(t). 
The convergence is found to be slower for a non-uniform 
initial tracer seeding. 

In the pair dispersion analysis, we confirmed the initial 
ballistic regime (R 2 ~ t 2 ) and a very short-range regime 
with R 2 ~ t 3 which have also been found in other studies 
on convection [TEl . [iH ] and simulations of homogeneous 
isotropic turbulence in a periodic cube (2(| [2{| . The fi- 
nite size of the convection cell suppresses the dispersion 
limit; hence a Taylor-like regime (R 2 ~ t) is absent for 
T < 8 in our studies as opposed to [l5|. The autocor- 
relations of the vertical velocity component are strongly 
anti-correlated. This is in line with [151 ] and seems to be 
associated with the characteristic structures - the ther- 
mal plumes - present in convection. These structures 
have a strong influence on the vertical tracer motion. In 



particular they cause the reversals in the tracer tracks, 
namely when a particle hits the top or bottom plates. 
While we did not observe a Rayleigh-number-dependence 
of the acceleration statistics, a systematic dependence of 
the PDFs of the acceleration components on T was ob- 
tained. To conclude, the basic Lagrangian properties we 
studied in the present work agree qualitatively with those 
in previous studies, e.g. in [15l |16I|. The additional vari- 
ation of the aspect ratio is in line with changes in the 
large-scale circulation in the convection cell, which af- 
fects the Lagrangian tracer dispersion and acceleration. 

Our studies complemented the only existing La- 
grangian laboratory experiment on convection by Gas- 
teuil et al. [lij]. We could confirm qualitatively some 
results obtained in their experiments. Further experi- 
mental and numerical studies in the large aspect ratio 
regime are desirable and necessary. Similar to the Eule- 
rian case, we can expect the aspect ratio dependence to 
disappear for sufficiently large values of F. 
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